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Abstract 

The nonthermal particle content of galaxy clusters should in part have a cosmolog- 
ical component generated during the early starburst phase of the member galaxies. 
This is reviewed in the framework of a simple cluster formation model suggested 
previously. It implies a nonthermal energy fraction of about 10 percent for the Intr- 
acluster gas. We also propose a mechanism for the early generation of Intracluster 
magnetic fields in terms of Galactic Winds. It results in typical field strengths of 
Gauss. Such comparatively weak fields are consistent with an inverse Comp- 
ton origin of the excess EUV and hard X-ray emission of the Coma cluster, given 
the radio synchrotron emission. The required relativistic electrons must have been 
accelerated rather recently, less than a few billion years ago, presumably in cluster 
accretion shocks. This is in contrast to the hadronic nonthermal component which 
accumulates on cosmological time scales, and whose 7r°-decay TeV 7-ray emission is 
expected to be larger, or of the same order as the inverse Compton TeV emission. 
This 7-radiation characterizes the energetic history of cluster formation and should 
be observable with future arrays of imaging atmospheric Cherenkov telescopes. 



1 Introduction 



Clusters of galaxies are the largest gravitationally bound structures in the Uni- 
verse and may confine a representative fraction of its mass. Therefore the study 
of their dynamical properties and radiation content should allow, amongst 
other things, conclusions on the relative amounts of baryonic and nonbaryonic 
matter in cosmology (e.g. White and Fabian, 1995, and references therein). 

Another basic characteristic, also due to confinement, is the ratio of thermal to 
nonthermal energy in these objects. To a significant extent that ratio is estab- 
lished during the epoch of galaxy formation and thus preserves the energetic 
history of cluster formation. We shall review this topic here. The confinement 
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Fig. 1. Large scale contour plot of an X-ray image of the Perseus cluster from the 
ROSAT all sky survey (courtesy H. Bohringer). 

of nonthermal particle components is intimately related to the existence of 
strong and chaotic magnetic fields in the intracluster medium (ICM), and we 
shall propose a mechanism for their early generation. This is followed by a 
discussion of the present-day nonthermal radiation from clusters in various 
wavelength ranges, in particular at very high 7-ray energies. 

Rich Clusters 

Rich clusters, i.e. conglomerates with typically more than 100 member galax- 
ies, have typical radii R c \ ~ few Mpc and baryonic masses M c \ ~ 10 14 to 10 15 M 
The closest large cluster is the Virgo cluster at a distance of d ~ 20 Mpc; all the 
others are at a distance d of hundred Mpc and beyond. Examples for bright 
and relatively nearby clusters are the Perseus and the Coma clusters with 
d ~ 100 Mpc. The Perseus cluster is the brightest cluster in soft X-rays (see 
Fig. 1). The large X-ray luminosity is due to the very hot (T ~ 10 7 to 10 8 K), 
massive (M gas ~ few x X;M ga i), and metal-rich ([Fe] c i ~ O.35[Fe] ) ICM gas 
(e.g. Bohringer, 1996). 

Apart from their primary cosmo logical interest, galaxy clusters also serve as 
extragalactic distance poles (Sunyaev-Zeldovich effect) and as gravitational 
telescopes for still more distant objects. This makes them important "instru- 
ments" for observational astronomy. 

Cluster Formation 

Most rich clusters - apart from objects still in an early formation phase, like 
the Virgo cluster - are characterized by a predominance of early type galaxies, 
i.e. elliptical (E) and lenticular (SO) galaxies which have little interstellar gas. 
Gas-rich spiral and irregular galaxies represent a minority. This situation is 
opposite to the one in the field, the space outside clusters, where the fraction 
of spiral galaxies is at least 75 percent. For rich clusters the metallicity of the 
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ICM gas, for instance in the form of the fractional ICM iron mass, is also 
correlated with the optical luminosity in E and SO galaxies (Arnaud et al., 
1992). The correlation supports the qualitative view that early starbursts due 
to galaxy-galaxy interactions of protospirals have produced a large number 
of Supernovae (SNe) that heated the originally present interstellar gas and 
generated violent Galactic Winds, leaving gas-poor E and SO galaxies behind. 
This mass loss should have led ultimately to the observed strong chemical 
enrichment of the ICM gas. We also conjecture that the radiation and the 
winds from these early galaxy mergers strongly heated the remaining primor- 
dial ICM gas and thus prevented further galaxy formation. This is perhaps the 
physical explanation for the observed inefficiency of galaxy formation which 
manifests itself in the remarkable preponderance of diffuse ICM gas mass over 
stellar mass in clusters of galaxies. 

A quantitative discussion of the dynamical prerequisites for Galactic Winds 
and the total number of SNe in clusters is contained in the paper by Volk et 
al. (1996, hereafter referred to as Paper I) which we shall summarize below. 

The total number of SNe since galaxy formation in the cluster, roughly a 
Hubble time T H ~ 1.5 x 10 10 yr ago, is given by 



iV SN = J dtx u SN (t) 



0.35 [Fe] x M cl 





5M Fe 

-T H 



where SM Fe is the amount of iron produced per event. In such starbursts 
we dominantly expect core collapse SNII from massive progenitor stars to 
occur, with 5M Fe ~ 0.1M Q on average. For the Perseus cluster this implies 
^SN rseuS ~ 3 x 10 12 . The corresponding total energy input into the interstellar 
medium is Asn-^sn ~ 3 x 10 63 i?5i erg, where E 5 i = 10 51 erg is the average 
hydrodynamic energy release per SN. 



Assuming the early starbursts to occur at a typical redshift of z ~ 2 due to the 
merging of protospirals in the overdense protocluster environment (Steinmetz, 
1993), with a duration of Tge < 10 9 yr, we obtain 

/ A /"Perseus / A /"Perseus \ 

^ ^ 1 > 100 x ^ lkyWay . 

Tsb 

Here N^{ seus ~ 500 denotes the number of galaxies in the cluster. As an 
example we can compare to the archetypical contemporary starburst galaxy 
M82. It has a current SN rate u^ 82 ~ 10 x z/^ ilkyWay , a wind velocity v wind ~ 
2300 km/sec, and a mass-loss rate of M ~ 0.8M Q /yr (Breitschwerdt, 1994). 
The starburst nucleus of M82 is characterized by the following values for 
the interstellar gas temperature T, gas density n, and thermal gas pressure p : 
Tbase ~ 10 8 K, n base - 0.3 cm" 3 , and p gas /k B ~ 10 7 Kcm~ 3 (Schaaf et al. 1989). 
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Since the thermal ICM gas pressure in the Perseus cluster is p^ rscus /kB ~ 
10 4 K cm -3 , it is clear that an object like M82 could readily drive a wind even 
against the present-day ICM pressure. At the galaxy formation epoch the ICM 
pressure was certainly much smaller than this value. 

In the expanding wind flow the SN-heated gas cools adiabatically to quite 
small temperatures. However it is reheated in the termination shock, where 
the ram pressure of the wind adjusts to the ICM pressure. Beyond this point 
the ejected galactic gas is mixed with the unprocessed ICM gas. 

Particle Acceleration 

Cluster formation also implies the production of a strong nonthermal compo- 
nent of relativistic particles. They will be accelerated during the early phase - 
and possibly also in later events - and confined in the turbulent ICM medium. 
The confinement time generally exceeds the cluster lifetime. Thus the energy 
spectrum of the energetic particles is the same as that generated by their 
sources. In Cosmic Ray parlance these are cosmological Cosmic Rays (CRs). 

During the early starburst particle acceleration will occur initially at the outer 
shocks of the Supernova Remnants (SNRs). However, like the thermal gas, this 
first generation of nonthermal particles will loose its energy almost completely 
by adiabatic cooling in the ensuing Galactic Wind. 

Fresh particle acceleration will occur with high efficiency at the strong wind 
termination shock, at distances ~ 100 kpc. The wind magnetic field will still 
intersect the shock at an angle ~10° (e.g. Zirakashvili et al. 1996), so that the 
standard process of diffusive shock acceleration approximately still works there 
(Fig. 2a). Assuming an overall acceleration efficiency of 10 to 30 percent, one 
gets a total gas internal energy ~ few 10 62 erg and a nonthermal energy 
~ 10 62 erg for a system like the Perseus cluster, ultimately from star for- 
mation and subsequent SN explosions. Since the galaxies are distributed across 
the cluster quasi-uniformly, this will also be true for the nonthermal particle 
population and the radiation they emit. The continuing gravitational contrac- 
tion/accretion of the cluster will subsequently energize CRs and thermal gas at 
least adiabatically, or shock accelerate/heat both components, so that finally 
the total energy E CR of energetic particles reaches 10 63 erg~ 0.1-E gas in the 
cluster; E gas now denotes the total internal energy of the ICM gas (Paper I). 
The resulting nonthermal energy density of some tenths of eV/cm 3 happens 
to be roughly equal to that in the interstellar medium of our Galaxy. 

It is instructive to compare the expected nonthermal energy with the thermal 
energy content of the cluster galaxies. Assuming the stars internally to be in 
virial equilibrium and, for purposes of estimate, all of them to have a solar 
mass and radius, then E^ ~ (3/lO)GM|/_R ~ 10 48 erg. For a total mass 
of about 1O 14 M contained in the galaxies of the Perseus cluster this gives 
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Fig. 2. a (left panel). Schematic of particle acceleration at a Galactic Wind termi- 
nation shock. In the halo, at distances large compared to the radius of the galaxy 
(Disk), the wind velocity becomes radial and constant in magnitude. It goes through 
a strong shock transition at a distance where the ram pressure of the wind has 
decreased towards the external (IC) pressure. Particles injected from the heated 
thermal gas of the shocked wind get diffusively accelerated at this shock with high 
efficiency b (right panel) . Schematic of the flow and field configuration of an axisym- 
metric Galactic Wind from a rotating disk galaxy like the Milky Way, cf. Ptuskin 
et al. (1997). The meridional flow velocity u from a circle in the disk (like that 
corresponding to the Sun's distance from the center of rotation) is roughly parallel 
to the axis of rotation to flare out at large distances. The magnetic field B becomes 
slowly azimuthal at large meridional distances s. 

a total thermal energy in stars ~ 10 62 erg, and thus E CR > X) ga i Sstars ^th ar - 
This means that the nonthermal ICM energy is larger than the total thermal 
energy of all the stars in all the galaxies contained in the cluster! 

It has been argued more recently that, apart from star formation and overall 
gravitational contraction, also individual giant radio galaxies should have in- 
jected large and in fact comparable amounts of nonthermal particles during 
the life time of a cluster (EnBlin et al. 1997; Berezinsky et al. 1997). This is 
no doubt an important additional possibility. A weakness of this argument 
consists in the fact that per se it is predicated on statistical knowledge about 
the luminosity function for active galaxies in clusters in general, and not on 
direct observations of the individual cluster to which it is applied. 

Particle Confinement 

The large-scale magnetic field in the ICM gas may be quite chaotic and not 
well connected over distances exceeding typical intergalactic distances (see 
section 2). Thus energetic particles may not readily escape from the cluster 
due to such topological characteristics. However, already pure pitch angle dif- 
fusion along magnetic field lines with superposed turbulent fluctuations gives 
important insights into the confinement properties of galaxy clusters. Stan- 
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dard quasilinear theory yields a spatial diffusion coefficient k» along the large 
scale field B due to a power spectrum P(k) of magnetic field fluctuations with 
wavelength A = 2n/k as 

B 2 

-II = (V3)cr g (p) /rdA;/p(A;r 

where p denotes particle momentum, kr g (p) ~ 1, and k denotes wavenumber 
of the field fluctuations. Let us assume a relative fluctuation field strength 
of order unity at the inter-galaxy distance 1/fco, i-e. a totally turbulent field 
P(k ) x ko ~ B 2 on this scale, and a power law form of P(k) = P(k )(k/k )~ n . 
Then the diffusion time across the cluster T esc ~ k\\ > T H for (cp) protons < 
10 17 eV and < 10 15 eV, for n = 3/2 and n = 5/3, respectively (Paper I). Also 
tf r s ° s tons 3> T H for nuclear collisions in the ICM gas. Therefore (except at sub- 
relativistic energies with their prevailing Coulomb losses), up to these energies 
CR hadrons accumulate in the cluster since the galaxy formation epoch, and 
that is what we called cosmological CRs before. The situation is different for 
relativistic electrons, which suffer radiative losses: energetic electrons observed 
now must either be secondaries or be rather recently accelerated. 



2 Intracluster Magnetic Fields 

The large magnetic field strengths of B ~ 1 /xG in the IC medium of rich 
clusters, in particular as observed by Faraday rotation measurements (e.g. 
Kronberg 1994), are not easily explained by a contemporary mechanism be- 
cause present day turbulent dynamo effects in such a large-scale system should 
be extremely slow. Therefore we suggest here a field configuration that is due 
to the early formation history of galaxy clusters and that should be essentially 
preserved to this day. It should even be still in a state of development at the 
present epoch. The argument derives from the violent early Galactic Winds 
which accompany the starbursts responsible for the predominance of the early 
type galaxies in rich clusters. 

We assume first of all that the protospirals, whose mergers constitute the build- 
ing blocks for the E and SO galaxies, had already generated galactic magnetic 
fields of /iG strength. This should indeed be possible within about 10 8 yr, i.e. 
of the order of a rotation period of our Galaxy, from turbulent dynamo action 
that invokes boyancy effects from CRs and magnetic reconnection on spatial 
scales of O(lOOpc) (Parker, 1992); such a time scale and the resulting field 
strengths correspond to generally accepted numbers. In the second stage the 
ensuing Galactic Winds extend these fields from the interacting galaxies to 
almost intergalactic distances. In the final and by far longest stage, that lasts 
until now, the fields are recompressed by the contraction of the cluster to its 
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present size. 



The ICM fields do not reconnect on the intergalactic scale in a Hubble time. 
Consequently there is no need for a continuous regeneration of these fields 
since their formation. However, this also implies that a topologically connected 
overall ICM field will on average not be formed either, and that the ICM field 
is chaotic on a scale smaller or equal to the present intergalactic distance. 

In detail we draw on arguments we have in the past used for the field con- 
figuration in a Galactic Wind from our own Galaxy (Zirakashvili et al. 1996; 
Ptuskin et al. 1997). They are based on estimates of the relative amount of 
field line reconnection vs. the extension of galactic field lines by a wind to 
"infinity" (Breitschwerdt et al. 1993). The basic result was that the rates 
of reconnection - and thus of the formation of "Parker bubbles" leaving the 
galaxy by their boyancy and allowing the generation of the disk magnetic field 
- and of extension of this field into the galactic Halo by the pressure forces of 
the wind are roughly equal. Thus both effects occur, and in the cluster galax- 
ies, magnetic energy can be generated on the large scale of the wind at the 
expense of the thermal and nonthermal enthalpies produced in the starburst. 
The geometry of the field should roughly correspond to straight field lines out 
to radial distances r of the order of the starburst (SB) radius, ~ 1 kpc, and 
spherically diverging field lines beyond that. The slow rotation of the system 
should then lead to an azimuthal field component oc 1/r which dominates at 
large distances over any radial component. However, in contrast to the famil- 
iar situation in the Solar Wind equatorial plane, the axis of rotation is rather 
parallel than perpendicular to the flow at the base of the wind, and thus the 
dominance of the azimuthal field component is by no means as drastic as in 
the case of a stellar wind (Fig. 3). 

Assuming the ICM pressure p c i(z = 2) at the formation stage of the early 
type galaxies to be roughly a factor of 10 -2 smaller than it is at present (after 
gravitational compression of the ICM gas), the termination shock distance r s h 
is found from p{r sil )u 2 /kB ~ 10~ 2 p c i(z = 0)/k B ~ 10 2 K/cm~ 3 . 

For M82 analogs, but with ~ 1 kpc, we have p gas /kB ~ 10 7 Kcm~ 3 . Thus 
~ Pgas/[10~ 2 Pci(0)] ~ 300, and therefore 

r sh 300 kpc 

~ U.o 



d^(0)/(l + z) (2Mpc/3) 

The Wind Bubble containing the shock-heated wind gas will have a radius still 
exceeding r s h. Thus, even though the volume of hydrodynamically unaffected 
ICM gas may be large enough so that the ICM gas mass exceeds the mass 
associated with galaxies by a factor of a few - as observed - there may be 
that some Wind Bubbles touch. However, at the scale of r sh , reconnection 
with a speed between 1 and 10 percent of the Alfven velocity is too slow to 
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occur over a Hubble time, even in a present-day ICM magnetic field as high as 
10~ 6 G. Therefore, on average, the field structure should remain topologically 
disconnected until today. 

The field strength B d (z = 2) in the Wind Bubbles should be of the order of 

B cl (z = 2) ~ 4S gal r^/r sh ~ 10- 2 5 gal ~ 1(T 8 G 

or somewhat larger, if the field in the bubble increases in the decelerating 
postshock flow. The ongoing cluster contraction/accretion compresses the field 
to lowest order isotropically oc I 2 , with the scale factor 

Z^[n cl (0)/n bar (0)] 1 / 3 /(l + ^), 

where n c \(0) ~ (10 _3 to 10~ 4 ) cm" 3 and z — 2. Choosing for the present mean 
baryon number density the value nb ar ~ 3 x 1CT 7 cm" 3 , we obtain / ~ (2.3 to 5). 

Finally then, we obtain for the present-day ICM field: B c \(z = 0)/B c \(z = 
2) ~ (5 to 25). Therefore the present-day ICM magnetic field should have a 
mean strength of the order of 10 _7 G, from "primordial" seed fields, randomly 
directed on an intergalactic scale. Although smaller by about one order of mag- 
nitude than estimated from Farady rotation measurements, such fields need 
not necessarily be unrealistic, considering that observations might emphasize 
regions of high magnetic fields. In addition the increase of the field in the 
Galactic Wind bubbles beyond their postshock value might be more than a 
factor of unity as assumed above. Thus we cannot exclude /iG fields although 
they certainly are at the upper limit our estimate permits. The interpretation 
of recent UV data also points to small field values (see section 3). 

In conclusion, there is hardly any need for a contemporary "turbulent IC 
dynamo". Nevertheless, the estimated CR enthalpy of ~ 0.5erg/cm 3 is essen- 
tially a free energy reservoir for a future increase in B towards a mean strength 
of a few fiG. In the next section we consider the nonthermal radiation from 
clusters. In particular we shall discuss the implications of the broad-band ob- 
servations for the nonthermal energy content and ICM magnetic fields. 



3 Nonthermal emission 

From the Coma cluster radio fluxes are measured from 10.3 MHz to 2.7 GHz 
(Bridle and Purton, 1968; Henning 1989; Kim et al. 1990; Giovannini et al. 
1993), and at 4.85 GHz an upper flux limit was reported by Kim et al. (1990). 
The energy fluxes J v oc z/~ Qr between 30.9 MHz and 1.4 GHz are well fitted 
with a power-law index a r = 1.16 (Bowyer & Berghofer 1998). The data at 2.7 
and 4.85 GHz fall below this extrapolation, which may indicate a steepening in 
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the electron spectrum (Schlickeiser et al. 1987) for Lorentz-factors 7 > 5 x 10 4 , 
but can be explained also as an instrumental effect (Deiss 1997). Below we 
use an index a r = 1.16. 

In the extreme ultraviolet (EUV) region, diffuse radiation between 65 and 
245 eV is observed which appears in excess of the thermal fluxes of the X-ray 
emitting gas with T ~ 2 x 10 6 K, and this radiation was initially interpreted in 
terms of thermal emission from a gas of lower temperature, T ~ 8 x 10 5 K (Lieu 
et al. 1996). Subsequently, the excess EUV emission, observed in a number of 
clusters (e.g. Mittaz et al. 1998), was suggested to represent inverse Compton 
(IC) radiation (Hwang 1997; Sarazin and Lieu 1998) of low energy electrons, 
7 ~ 300, on the 2.7 K microwave background radiation. Also EnBlin and 
Biermann (1998) used this possibility to estimate lower limits to the magnetic 
field strength in the Coma cluster (see also Lieu et al. 1998). 

Bowyer and Berghofer (1998) have shown that the size of diffuse emission 
of the Coma cluster is significantly larger in radio than in EUV light. In 
addition, the spectral index of the energy flux J(E) of the EUV emission 
is a uv pa 0.75 which is significantly different from that in the radio domain. 
Based on these arguments, Bowyer and Berghofer (1998) suggested that two 
different populations of relativistic electrons should be responsible for the radio 
synchrotron and the EUV IC fluxes. We shall show below that the differences 
in these spectral indices are quite naturally explained in terms of a single 
population of electrons producing both the EUV excess IC and the radio 
synchrotron emission which avoids the need for a 2-component model. 

The characteristic synchrotron frequency produced by electrons with Lorentz 
factor 7 in the magnetic field B is v ~ (B/l /iG)7 2 Hz (e.g. Ginzburg 1979). 
Thus in a magnetic field of B ~ 10 -7 G the production of synchrotron ra- 
diation with a T = 1.16 in the region from 30 MHz to 1.4 GHz requires the 
energy distribution of the electrons ^(7) oc 7"°° to extend with an index 
tt e = 1 + 2a r pa 3.3 from 71 ~ 1.7 x 10 4 to 72 ~ 1.4 x 10 5 . These Lorentz 
factors are smaller by a factor of 3 if B ~ 10~ 6 G. 

For the IC process the mean energy of photons is e ~ (4/3)e 7 2 , where e 
is the energy of target photons. For the 2.7 K background radiation eo = 
6.5 x 10~ 4 eV. For the production of Coma's EUV spectrum from 65 to 250 
eV one therefore needs a spectral index of electrons a c ~ 2.5 only in a rather 
narrow energy region around 7 ic ~ 300. This energy is much smaller than 
the energy 7x necessary for the radio emitting electrons. Thus a scenario in 
which radiative losses induce a break in a e for a single population of electrons 
somewhere in between 7; c and 71 can readily account for the observed indices 
both in the radio and the EUV. Indeed, assuming that the spectral index of 
the source function of electrons in the ICM is = 2.3, the IC emission of 
those electrons at low energies will give a uv = 0.65, which is rather close to 
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Fig. 3. Synchrotron and IC fluxes (top panel) and the energy distribution of elec- 
trons (bottom panel) calculated assuming continuous (heavy lines) and impulsive 
(thin lines) injection of electrons with a- m j = 2.3 and j c = 2 x 10 8 , for the in- 
jection times t = 10 9 yr (solid) and t = 10 10 yr (dashed). The ICM gas density 
n = 10~ 3 cm~ 3 . For the assumed B = 0.1 //G the electrons are stationarily injected 
with Linj = 2.6 x 10 45 erg/s; the total energy input during 1 Gyr is 8.2 x 10 62 erg. 
For the impulsive injection, the same energy inputs as in the corresponding cases of 
continuous injection are assumed. Nonthermal fluxes from Coma in the radio, EUV, 
X-ray, and the EGRET upper flux limit (Sreekumar et al. 1996) are also shown. 

the value 0.75, in a rather narrow EUV band. After the break, the electron 
distribution steepens to a e = 1 + a in j = 3.3, required for the radio emission of 
Coma. 

Such a break in the energy spectrum is necessarily produced by radiative losses 
of electrons on time scales > 10 9 yr. Indeed, the characteristic energy loss time 
in the 2.7 K background radiation field as well as due to synchrotron emission 
can be written as 

t rad = 2.4 x 10 12 [C( 7 ) + 0.LB 2 G rV 1 yr , (1) 

where B^ G = B/l jj,G. The coefficient Cij) takes into account corrections for 
the Klein-Nishina effect in the IC energy losses of ultrarelativistic electrons, 
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which is important for B < 3 fiG. C{^) = 1 for 7 <C 10 s (IC losses in the 
Thompson limit), but C(10 8 ) ~ 1/2, C(2 x 10 8 ) ~ 1/3, and C(5 x 10 8 ) « 1/6. 

For a magnetic field < 3 /iG and timescales ~ (1 — 3) x 10 9 yr, equation (1) 
predicts the radiative cutoff energy 7 br ~ 10 3 . It is important to note that 
relativistic electrons need to be produced (accelerated) in the ICM continu- 
ously during all these last years, because in the case of an 'impulsive' injection 
of the electrons the radiative losses would remove all particles with energies 
above 7b r , and then the radio spectra cannot be explained. 

This is seen in Fig. 3 which shows the results of calculations that assume both 
continuous and impulsive injection of electrons with a spectrum 

Q(7) 0C7^exp(- 7 / 7c ), (2) 

where 7 C defines the assumed characteristic maximum energies of accelerated 
particles. The fluxes for continuous injection are normalized to the radio flux 
2 Jy observed at 400 MHz. The IC radiation for continuous injection during 
the last t = 10 9 yrs (heavy solid line) has the spectral shape of the EUV flux, 
and for a slightly smaller magnetic field (cf. Fig. 4), or an injection time larger 
by a factor of 2, it can also explain the absolute EUV flux. However in the 
case of relativistic electron production on much larger time scales, ~ 10 10 yr, 
the explanation of the shape of the EUV radiation becomes problematic even 
for continuos injection. Moreover, the agreement with observed spectra be- 
comes worse than shown in Fig. 3 if we take into account that for cosmological 
timescales the energy density of the microwave background increases with red- 
shift T A (z) oc (1 + z) A . Thus, the electrons responsible for the excess 
EUV radiation, if indeed it has an IC origin, must have been produced con- 
tinuously during the recent (1-3) Gyrs. This would be the characteristic age 
of the accretion shocks in the cluster which seem to be the most probable 
accelerators for the radio emitting electrons. 

An IC origin of the excess EUV flux imposes a strong constraint on the mag- 
netic field in Coma. The variation of the IC flux for different magnetic fields 
are shown in Fig. 4 which demonstrates that the lowest consistent magnetic 
field strength in Coma is equal to B — 7.5 x 10~ 8 G. This field could explain 
also the excess flux observed by Beppo-SAX in hard X-rays beyond 25 keV 
(Fusco-Femiano et al. 1998), that might be due to IC radiation of GeV elec- 
trons. The dot-dashed curve in Fig. 4 falls below both the X-ray and the EUV 
data. However, the assumption of injection for about 5 Gyr shifts the position 
of the radiative cutoff energy by factor of 5 and increases by the same factor 
the number of low energy electrons. The number of electrons with 7 ~ 300 is 
also increased for a steeper injection spectrum, a- m j = 2.5. Both these options 
could marginally explain the radio synchrotron and EUV IC fluxes. Thus an 
IC origin of the excess EUV radiation restricts the cluster magnetic field to a 
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Fig. 4. Nonthermal radiation fluxes calculated for continuous injection of electrons 
during t = 10 9 yr, for different ICM magnetic fields: B = 7.5 x 1CT 8 G (solid lines), 
1.5 x 1(T 7 G (dot-dashed lines), and 3 x 1(T 7 G (3-dot-dashed lines) 

narrow range around 0.1 /iG. 

For magnetic fields B ~ 10 _7 G the observed radio fluxes require injection of 
accelerated particles with a luminosity L in j ~ 3 x 10 45 erg/s. In a thermal gas 
with density n ~ 10~ 3 cm -3 , a secondary origin of the relativistic electrons, i.e. 
their production in pp interactions, would require an enormously large total 
energy in relativistic protons. Indeed, the luminosity in the — /i ± -decay 
electrons can be estimated as 

L± ~ 7.7 x 10 40 (n p /l(T 3 cm- 3 )- 3 (£ CR /10 60 erg) erg/s. (3) 

The fluxes of IC gamma-rays to be expected at TeV energies are shown in 
Fig. 5. for the case of B = 1CT 7 G and acceleration of electrons up to an 
exponential cutoff energy E c = 100 TeV. For magnetic fields ~ 0.1 /xG the 
diffusive shock acceleration mechanism in the Bohm limit still allows such 
high values for accelerated electrons. In this limit the acceleration time, defined 
by the rate d7/dt ~ ^u 2 / D(^) where u is the shock speed and ^(7) is the 
diffusion coefficient, is estimated as 

^ c ~5.4xlO- 5 % -lB^7yr, (4) 

where u 3 = -u/10 3 km/s. Equating t acc with t rad shows that values 7 > 10 8 are 
possible. 

The solid line in Fig. 5 corresponds to the solid curve in Fig. 3. The TeV fluxes 
of IC 7-rays in that case correspond to several per cent of the TeV emission 
from the Crab Nebula. In the case of single power law injection with ctinj = 2.3 
a further increase of the IC fluxes at TeV energies are impossible because 
of the EGRET upper limit above 100 MeV (Sreekumar et al 1996). Note, 
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Fig. 5. IC 7-ray fluxes expected from the Coma for B = 7.5 x 10~ 8 G and = 2.3 
(solid lines), and B = 9 x 10~ 8 G and a- m j = 2.1 (dashed lines). The dots show the 
7-ray fluxes produced in pp interactions of CRs with Eqr = 3 x 10 62 erg in ICM 
with n = 10 -3 cm~ 3 . 

however, that for a more realistic modelling of Coma as a spatially nonuniform 
source, the power law index a in j could be somewhat smaller than 2.3. We do 
not consider here this possibility, which should address also the question of 
different angular sizes of Coma in EUV and radio, as well as the flattening 
of the radio spectra observed towards the core of the cluster (Giovannini et 
al. 1993; Deiss et al. 1997). In Fig. 5 the dashed line shows the maximum of 
the IC radiation fluxes to be expected in the case of a in j = 2.1, which is then 
at the 0.1 Crab level. The dots correspond to the 7r°-decay fluxes produced 
by hadronic CRs with a cr = 2.1, assuming Eqr = 3 x 10 62 erg. The size of 
the radio emission produced by high-energy electrons in Coma is about half 
a degree. Therefore we can expect in that case a similarly large size for the 
TeV emission. Detection of such extended and weak (<0.1 Crab) TeV fluxes 
by current instruments is problematic, but they are quite accessible for the 
future HESS and VERITAS arrays. 

These results for the IC fluxes relate to the case of ICM magnetic fields 
B ~ 10~ 7 G. However, the estimates for the magnetic fields in Coma range 
from 0.1 to several /xG. In particular, Kim et al. (1990) and Feretti et al. 
(1995) deduced magnetic fields B ~ 1.7/xG and B ~ 6.0 fiG, respectively, 
from Faraday rotation measurements on background radio sources. If this is 
so, then an IC origin of the excess EUV fluxes is absolutely excluded. The only 
other possibility left for an explanation of this radiation in terms of nonther- 
mal radiation is synchrotron production by very high energy electrons. The 
interpretation of the steep radio fluxes then requires a second component of 
nonthermal electrons. 

For a production of synchrotron photons with energies 200 eV [y ~ 5 x 10 16 Hz) 
in magnetic fields B ~ 1 fiG one needs electrons with 7 > 2 x 10 s . From 
Eq.(l), the radiative loss times of these electrons is very short, t ra< j < 10 4 yr. 
Even assuming a rectilinear propagation of such electrons, the maximum pos- 
sible distance from the acceleration sites would be less than 3kpc. Thus, the 
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acceleration sites of the electrons should be rather smoothly distributed in the 
ICM in order to result in a smooth distribution of the synchrotron EUV ra- 
diation. This seems quite reasonable for accretion shocks, but not for isolated 
point sources like active galaxies. Another requirement is that the spectrum 
of accelerated particles should be very hard, with a in j ~ 1.5. Then the spec- 
trum of electrons steepens to a c ~ 2.5 which is required for the production of 
synchrotron radiation with a uv ~ 0.75. 

In Fig. 6 we show the spectra of the synchrotron and IC radiation generated 
in a magnetic field B = 2 jiG by two different populations of relativistic elec- 
trons. The dashed curves are produced, as in previous figures, by electrons 
continuously injected into the ICM during recent epochs up to t > 10 9 yr with 
a in j = 2.3 and luminosity L in j = 3.5 x 10 42 erg/s. For this luminosity, these 
electrons can be well explained as secondaries produced in the interactions of 
hadronic CRs, of total energy E C r = 4.5 x 10 61 erg/s, with an ICM gas of 
density n p = 10~ 3 cm~ 3 . They can generate the observed radio emission. On 
the other hand, the EUV fluxes in Figs. 3 and 4 are produced by the second 
component with 7 C = 6 x 10 8 , and a hard injection spectrum with a in j = 1 
(solid and dot-dashed curves) and a in j = 1.5 (three-dot-dashed). The curves 
show that continuous injection of this second component over times exceeding 
3 x 10 7 yr from now, without reacceleration, would lead to an accumulation of 
electrons at energies smaller than 7 ~ 10 5 . This would be quite sufficient to 
produce a radio synchrotron flux above the observed level. Moreover, in the 
case of a in j = 1.5 one has to assume a cutoff in the injection spectrum below 
7 low ~ 10 5 . Otherwise the number of low energy electrons in the injection 
spectrum would be unacceptably high. 

Thus, the scenario with high magnetic fields and nonthermal origin of the 
excess EUV radiation requires that the second electron component enters a 
reacceleration cycle at least once in 3 x 10 7 year. In addition, this compo- 
nent should not contain electrons with low energies. Therefore it cannot be 
produced by reacceleration of secondary electrons. 

In a speculative vain this component could be due to run-away pulsars born 
with kick velocities ~ 1000 km/s. Another possibility for the pulsars to appear 
in the ICM would be that they could be dragged there in the process of galaxy- 
galaxy collisions. Pulsars produce, at the pulsar wind termination shocks, an 
electron distribution with a strong deficit of low-energy particles (e.g. Arons 
1996). These electrons could be kept at those high energies by entering into 
frequent reacceleration processes on the accretion shocks in the ICM, and 
could be thus responsible for the second component. 

Reacceleration would not affect the radio electrons if it happened in the central 
region (see Bowyer and Berghofer 1998). The signature of the second compo- 
nent would be a rising spectrum of IC 7-rays. Fields below 1 jiG would even 
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Fig. 6. Nonthermal fluxes expected in the case of B = 2 /iG in a two-component 
model for relativistic electrons in Coma. For the first (radio) component (dashed 
lines), continuous injection during t = 10 9 yr is assumed, with a spectrum of the 
form of Eq.(2) with a- ul ^ = 2.3. For the second (EUV) component the parameters 
of the continuous injection, without reacceleration of electrons, read: t = 10 9 yr and 
ainj = 1 (solid lines), t = 3 x 10 7 yr and a- m j = 1 (dot-dashed), t = 3 x 10 7 yr 
and a; n j = 1.5 (three-dot-dashed). For the case of = 1.5 a rapid transition to 
a oc 7 2 behavior in the injection spectrum in the region 7 < 10 5 is assumed. The 
exponential cutoff energy is j c = 6 x 10 8 , and the gas density is n = 10~ 3 cm~ 3 . 
Injection rates for the first and second components are L m j = 3.5 x 10 42 erg/s and 
Lj n j ~ 5.5 x 10 44 erg/s, respectively. The total number of electrons in the second 
component is N c ~ 10 58 . The heavy bar corresponds to the "0.1 Crab" flux level. 

imply uncomfortably large TeV fluxes from the center of the Coma cluster. 
Thus the alternative between low and high magnetic fields in the Coma cluster 
would be given by the shape of the TeV 7-ray spectrum. 
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